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Abstract 

Integrable inhomogeneous versions of the models like NLS, Toda chain, Ablowitz-Ladik model 
etc., though well known at the classical level, have never been investigated for their possible quan- 
tum extensions. We propose a unifying scheme for constructing and solving such quantum integrable 
inhomogeneous models including a novel inhomogeneous sine-Gordon model, which avoid the dif- 
ficulty related to the customary non-isospectral flow by introducing the inhomogcncitics through 
some central elements of the underlying algebra. 
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1. Introduction 

Over more than past two decades active interest has been focused in the study of inhomogeneous 
integrable models. Various types of space and time inhomogeneities have been introduced successfully 
in the well known classical integrable models like the nonlinear Schrodinger equation (NLS) through 
linear, quadratic, cylindrical, radial etc functions [1, 2, 3, 4], the Toda chain (TC) with arbitrary 
inhomogeneity [5] and the Ablowitz-Ladik (AL) model with linear in space and arbitrary in time 
functions [6, 7, 8] etc., preserving their integrability. The associated Lax operators, soliton solutions, 
Painleve integrability criteria etc. have been extensively studied for such systems [9, 10]. Such 
inhomogeneities might induce intriguing effects like locally varying interactions, space-time dependence 
of soliton velocities, trapping of solitons in periodic movements etc. In the simplest example of linearly 
inhomogeneous NLS the solitons move with a uniform acceleration [1]. 

However, it is rather surprising that, investigations for the above inhomogeneous models have 
been carried out only at the classical level and no systematic effort has been made toward their 
quantization, except perhaps our own preliminary study [11]. A possible reason for this might be, 
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that the inhomogeneities in such models are believed to be generated through non-isospectral flows 
with space-time dependent spectral parameter A = X(x,t). This in general would lead to a dynamical 
classical r-matrix and consequently to a quantum i?-matrix with space-time dependence, taking us 
beyond the known formulation of quantum integrable systems. 

We overcome this difficulty for some of the nonisospectral flows and propose a unifying quan- 
tization scheme for the integrable inhomogeneous NLS, TC and AL models together with a novel 
inhomogeneous sine-Gordon model. The idea of such quantization is based on our earlier scheme [12] 
as well as on the observation that, certain non-isospectral problems may be looked from a bit different 
angle by considering spectral parameter A still to be a constant, while relegating the inhomogeneities 
to a set of central elements of an underlying algebra, which ensures quantum integrability. These 
inhomogeneity elements, which may be interpreted as external classical fields, commute with all other 
basic operators and appear only in the quantum Lax operators, but not in the quantum .R-matrix. 
Such fields might have arbitrary space-time dependence with suitable restrictions on their boundary 
and asymptotic conditions, imposed by the integrability. A positive outcome of this approach is that, 
the classical and the quantum R matrices for such integrable inhomogeneous models remain the same 
as in their homogeneous counterparts. 

Note that the algebraic Bethe ansatz solution of quantum integrable systems depends on the 
vacuum Lax operator as well as on the i?-matrix elements. Therefore, since in our inhomogeneous 
extensions, the i?-matrices are kept the same, the related result are affected only by the inhomogeneity 
in their Lax operators, which signifies the presence of impurities, defects, density fluctuations in the 
media, or the influence of variable external fields, depending on the physical situations [16]. 

We concentrate here on inhomogeneous quantum NLS, and TC models as well as AL and SG 
models, in their different forms. Note that homogeneous versions of all these models are quantizable 
[13, 14], where the first two sets of models are linked to the rational, while the last two to the 
trigonometric i?-matrix. 

The arrangement of this paper is as follows. We briefly review in sect. 2 various known forms of 
the classical inhomogeneous equations and in sect. 3 the basic structures of the quantum integrable 
systems. We present in sect. 4 the construction of quantum integrable extensions of inhomogeneous 
models in a unifying way and in sect. 5 show their exact and systematic solution through the Bethe 
ansatz method. Sect. 6 gives the concluding remarks. 
2. Classical integrable inhomogeneous models 

We briefly list different well known forms of the classically integrable inhomogeneous NLS, 
TC and AL models, where the details can be found in the cited references. 
2.1 Inhomogeneous NLS equations 

1. NLS equation with linear inhomogeneity (XNLS) 

An integrable NLS equation (NLSE) with linear x-dependent inhomogeneity was proposed 
in [1, 2] with non-isospectral flow A = a, A(t) being the spectral parameter. The inhomogeneous 
equation allows uniformly accelerated soliton solution as well as time dependent wave-number 
and frequency of the enveloping wave [1]. 

2. Cylindrical NLSE (CNLS) 

An integrable cylindrically symmetric NLSE having an explicit time-dependent coefficient 
was proposed at the classical level, with non-isospectral dependence of spectral parameter: 
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A(x,t) = ^ + f t [4]. 

3. NLSE with t-dependent coupling (TNLS) 

It is shown through Painleve analysis that the time-dependence of the nonlinear coupling 
in the integrable NLSE can be of the form F(t) = This model also has the same non- 

isospectral flow as the previous case. 

4. Radially symmetric NLSE (RNLS) and 

5. x-dependent nonlocal NLSE (NLNLS) are known to be classically integrable [2, 4]. 
In both these cases the non-isospectral condition was taken to be A = aA 2 . 

6. NLSE with more general inhomogeneity (FNLS) 

Integrable inhomogeneous classical NLSE with more general x-dependent coefficient was 
proposed as iip t + V'xx + 2(|^| 2 — F(x))ip = 0, where F(x) can be linear, quadratic or in more 
general form, depending on the type of non-isospectrality given by the space-time dependent 
spectral parameter satisfying A t = 2(A 2 ) X F X [3]. 

2.2 Inhomogeneous Ablowitz-Ladik model 

The AL model was discovered first as a discrete NLSE. The inhomogeneous AL model was 
studied at the classical level in a series of papers [6, 7, 8], which describes the system to be 
in an external time- dependent and linear in space potential and induces an intriguing effect of 
trapping the soliton and forcing it to a periodic movement. 

2.3 Inhomogeneous Toda chain 

Classically integrable inhomogeneous Toda chain with varied space-time dependence, in of 
the form [5] u t t(n) = g 1 ( y n)e u ^ n ~ 1 ^~ u ^ — gi(n+ i) e u (™)- u ( n + 1 ) _|_ g 2 (n)+ boundary terms, allows 
different choices for inhomogeneity functions g a (n,t),a = 1,2. 
3. Quantum integrable systems 

As is well known [13], the quantum integrability of a system is guaranteed by the quantum 
Yang-Baxter equation (YBE) 

R(X - rfLjiX) ® Lj(n) = (I ® Lj(fj,))(Lj(X) ® I)R(X - »), (1) 

with the quantum Lax operator Lj of the discretized model defined at all sites j = 1,2, ... ,N 
and a site-independent i?-matrix, together with the ultralocality condition: (I ® Lk(fj,))(Lj(X) ® 
/) = Lj(X) <S> Lk(n), j ^ k. The corresponding classical model with the related r-matrix: 
R(X) -> I + hr(X) + 0(h 2 ) would satisfy the classical YBE {Lj(X) ®, L k (fi)} = 5 jk [r(X - 
fi),Lj(X) <E> Lj(fj)] for the discretized Lax operator Lj(X), which goes to the corresponding 
field Lax operator U(x, X) at the continuum limit A — > 0: Lj(X) — > / + iAU(x, A). In such an 
approach therefore only the space-Lax operator Lj(X) or its continuum version U (x, X) together 
with the canonical relation play the central role, while the time-Lax operator V(x, X) becomes 
insignificant. Moreover in quantum problems, unlike classical equations of motion, the main 
emphasis is on solving the eigenvalue problem of the Hamiltonian: H\n >= E n \n >, together 
with that of other conserved operators Cj,j = 1,2,.... These conserved operators are generated 
by the transfer matrix r(A) = tr(Y[j Lj(X)) through expansions r(A) = J2C±jX ±: > or through 
similar expansions of logr(A). 

The algebraic Bethe ansatz, an exact method for quantum integrable systems solves therefore 
a general eigenvalue problem: r(A)|n >= A n \n >. As an essential condition for the quantum 
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integrability , which follows from the ultralocality and the quantum YBE (1) the conserved 
operators including the Hamiltonian must commute mutually [13]. 

It is important to note that, the YBE (1) not only ensures the quantum integrability but 
also plays a central role in the exact Bethe ansatz (BA) solution and for this it is essential that, 
the i?-matrix appearing in it must not depend on the space-time coordinates. Consequently, 
since the R(X — //)-matrix is a function of A — //, the spectral parameter can not dependent 
on space-time variables, even in inhomogeneous models. Therefore, for quantization of our 
inhomogeneous models we introduce inhomogeneity through Lax operator Lj, by keeping the 
spectral parameter constant and demanding that the changed Lax operator must satisfy the 
YBE (1) with the same i?-matrix as in the original homogeneous model. A similar argument 
holds for the classical YBE at the classical limit. Therefore, it should be clear that we can 
consider inhomogeneous extensions for only those models, which are quantum integrable at 
their homogeneous limit. The inhomogeneous versions of integrable NLSE, TC and AL models 
discussed above are known only at the classical level as a non-isospectral problem, whereas 
inhomogeneous SG model seems to not have been studied even at the classical level. Our main 
concern now would be to look into the possible quantum extensions of these inhomogeneous 
models, which apparently have never been undertaken. 
4. Unifying quantization for inhomogeneous models 

We formulate a quantization scheme for inhomogeneous integrable models of both rational 
and trigonometric type based on an unifying quantum algebra [12]. The inhomogeneities are 
introduced in these models through central elements of the underlying algebra in a systematic 
way. 

4.1 Rational class of inhomogeneous models: 

Recall that both the homogeneous NLS and TC models are quantum integrable and associ- 
ated with the well known rational i?(A)-matrix [14] given by its nontrivial elements as 

a(A) = R{{ = Rf 2 = A + 77, b(\) = = R 2 2 \ = \,c = R% = R\\ = r), (2) 

For constructing the intended inhomogeneous extensions of these models we start from a rational 
discrete general Lax operator [12] 

r / \ \ / C ?(A + s^) + c] s n \ , 

4 c°(\-sl)-c 2 )> (3) 

which was shown to yield all integrable models with 2x2 Lax operators, belonging to the 
rational class. Note that L-operator (3) depends on operators s with generalized spin algebra 

[4> S 'J = S nm (2m + sl + m~), [si s±] = ±<5 nm (4) (4) 

and a set of commuting operators c", a = 0,1, a = 1,2, forming central elements m + = 
c\c 2 , mr = c\c 2 + c\c\, which commute with all other operators. We note that, algebra (4) 
dictated by YBE (1) guarantees the quantum integrability of all models realized through the 
Lax operator (3) and since the algebraic relations are valid locally at all sites n — 1,2,..., the 
central elements can be taken in general to be time as well as lattice-site dependent variables: 
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c a,nW- We would see in the sequel that these inhomogeneous elements would in fact be crucial 
in introducing inhomogeneity in integrable quantum systems we intend to construct. 
I. Quantum inhomogeneous NLS model 

It is easy to see that when inhomogeneous elements are absent i.e. when c\ = c° = c\ = 
—c\ = 1, giving m + = l,m~ = 0, the generalized spin algebra (4) reduces simply to su(2) spin 
algebra and through bosonic realization of its generators by the Holstein-Primakov representa- 
tion (HPR) we can recover from (3) the Lax operator L l ™ ls of the exact lattice version of the 
standard quantum NLS, introduced by Korepin and Izergin [15]. Therefore for constructing the 
corresponding inhomogeneous quantum model our strategy would be to follow a similar path, 
but with nontrivial elements c% n (t), for which we first find a realization of (4) in the form of 
generalized HPR 

s 3 n = s n -N n A, s+ = f(N n )i/> n y/A, s~ = ^J(N n )y/K, 
with f{N n ) =m~+ m+(2s n - N n A), N n = ^J) n . (5) 

where [tp n , = j^5 nm , s n is the spin parameter taken also to be site-dependent and A is 
the lattice spacing. (5) clearly recovers the standard HPR at the homogeneous limit. There- 
fore in analogy with [15] we show that, the generalized HPR yields an exact lattice version of 
the inhomogeneous quantum NLS model, for which we demonstrate first that, the associated 
L-operator (3) has a consistent continuum limit, which recovers the Lax operator of the corre- 
sponding quantum field model at A — > 0. Making a particular reduction of the central elements 
as 

Ci=4 = g n (t), c\ = -c\ = f n (t), giving m + = g 2 n , m~ = 0, (6) 

where f n ,g n are space-time dependent arbitrary functions, we find that, at the high spin limit 
s n — > \9n l the generalized HPR (5) reduces (3) to 

T lnlsUnh)(\\ \ _3 r ( \\ T,a( ~ 9nAN n tpnf^{N n )y/%, \ 

where A n = g n \ + f n and f^(N n ) = (2 — g n A 2 N n )^. One can verify that the related discrete 
model is a quantum integrable system, since the Lax operator (7) associated with it together 
with the rational i?-matrix (2) exactly satisfy the quantum YBE (1). We observe further that, 
at the field limit A — > 0, when ip n (t) — > ip(x,t), f n — > f(x,t), g n — > g(x,t), the lattice Lax 
operator (7) reduces to: L^(^) = I -iAU nls( - inh \x, A) + 0(A 2 ), recovering the corresponding 
field Lax operator 

«MH(_4^ A ) (8) 

where A = g(x,t)X + f(x,t). We may check again that, (8) with quantum field operators 
[${x,t),ift{y,t)\ = S(x — y) and f(x,t), g(x,t) acting as classical background fields satisfies 
YBE (1) with the same i?-matrix (2), up to the first order in 0(A), which however is sufficient 
for the integrability of a quantum field model. Therefore the quantum Lax operator (8) clearly 
represents an inhomogeneous generalization of the quantum integrable NLS field model. 
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For the homogeneous lattice NLS model it has been established [15] that, at certain values of 
the spectral parameter A = z/j, where the quantum determinant qdet(L(A)) = ^tr(L n (X)a 2 L n (\+ 
l)cx 2 ) vanishes the Lax operator is realizable as a projector as well as an inverse projector. This 
important property was shown to be useful for proving the local character of the discrete Hamil- 
tonian, i.e. its interaction spreading over to only few nearest neighbors. For our inhomogeneous 
lattice NLS model we could show that qdet(L(A)) for the associated Lax operator (7) vanishes 
at Vi = — ^, = ((? n A — ui) and the L n -operator at these points can be represented indeed 
as projectors and inverse projectors as 

Llnls(inh) / \ n i T Inls(inh) / (n)\ ~ n / n \ 

.. >(is 1 )=a j (3 i and L tj \v\>) = aft i (9) 

with P ± = <5f 2 = AipnyKh, «i = h = -A?Piy/(h, a 2 = p 2 = f i0) (N n ), Qfi = j3i = 
/ (0) (A n -l). 

However we note that the point becomes site-dependent for this inhomogeneous model 
and therefore the projector-representation can not be achieved for all L n at a single value of 
A, which would lead inevitably to nonlocal interactions for its Hamiltonian in both classical as 
well as quantum cases. By choosing some particular forms for g n , e.g., taking them same for 
all odd n — 2k + 1, we can partially regulate the localization of the Hamiltonian, though such 
inhomogeneous models in general would lead to nonlocal interactions. 

A further technique of pairwise grouping of the Lax operators at different sites for achiev- 
ing locality for the Hamiltonian has been shown for the homogeneous quantum lattice NLS. 
Whether the same technique is applicable also to our corresponding inhomogeneous case needs 
more detailed analysis , which we leave for future study. 

At the classical limit, when field operators become classical functions one can show that, 
(7) satisfies classical YBE with the classical rational r(A) = ^(J + a • <r)-matrix and represents 
therefore the Lax operator of a new integrable discretization of inhomogeneous classical NLSE. 
Similarly, at the classical limit (8) must also satisfy the continuum version of the classical 
YBE: {U(x, A) ® j U(y, fi)} = 5(x - y) [r(A - fi), U(x, A) <g> I + I <g) U(x, /i)], with the same 
r-matrix and would correspond to the classically integrable inhomogeneous NLS field models. 
Notice that, if in (8) one considers A to be the spectral parameter instead of A entering r(A)- 
matrix, the problem would look like a customary non-isospectral flow with nontrivial A t ,A x . 
In fact if we rename = Q, (8) would coincide at its classical limit with the Lax operator 

of known inhomogeneous NLS models and hence would recover for different choices of the 
functions f(x,t),g(x,t), the inhomogeneous NLS equations proposed earlier. For example, i) 
g — 1, f — at would give linear time- dependence of A = A + at with A = a reproducing 
the XNLS, while ii) g — |, f = f t would give the inhomogeneous CNLS. We can get also 
the inhomogeneous TNLS by multiplying the field variable in CNLS by a function of t, which 
however is not a canonical step. We may consider now a more general situation iii) g=X(x), 
f=T(t)X(x), which gives A = X(x)(X + T(t)) = g(x)h(t), coinciding with the separable solution 
of the FNLS [3]. Note that, the related non-isospectral picture satisfying the equation T = 
X'(T — A) 2 + F'(x)/X, is compatible with our isospectral relations. Therefore a nontrivial 
solution [3] : X = ax + b,T = a(T — A) 2 + a , yielding the NLS with quadratic x-dependence: 
F(x) = ao(|ax 2 + bx) + c, should also be reachable from our construction. 
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In cases of radial and nonlocal NLS however it appears that, due to nonlinear nature of the 
non-isospectrality and more involved noncanonical fields, the known classical models RNLS and 
NLNLS would not coincide with the classical limits of our quantum inhomogeneous models. 
It is important to note that, in classical Lax pair approach non-isospectrality with space-time 
dependence appears in both the space (U) and time (V) Lax-operators and moreover the allowed 
transformations may not respect the canonicity of the fields. In the quantum case on the other 
hand the inhomogeneity can appear only in the space part U or Lj, since the analog of V- 
operator is absent here. At the same time, the canonical structure defined in quantum models 
is fixed by the associated i?-matrix and should remain unchanged under all transformations. 
Nevertheless, as we have seen, apart from a simple noncanonical transformation of the field: 
tp — > Q = \^2gip, with g a space-time dependent function, (8) in the classical limit maps into 
the known Lax operator for the inhomogeneous NLS equations. 
II. Quantum inhomogeneous Toda chain 

We intend to construct the quantum extension of the inhomogeneous Toda chain start- 
ing again from the same Lax operator (3), but choosing now the inhomogeneity fields c\^c\j 
arbitrary, while c\ = c° = 0. This clearly gives both mf = reducing algebra (4) simply to 

[s+,s-] = 0, [sls±] = ±6 nm s± (10) 

which can be realized through canonical quantum fields [u n ,p m ] = i5 nm of the Toda chain as 



s 3 



) = -i Pj , sf = (G j f 1 e^, (11) 



with arbitrary function Gj. This reduces (3) to the quantum Lax operator of the inhomogeneous 
Toda chain: L T 3 C{inh \\) = L 3 (\)/c% as 

L > {X) -\ e^F m J' (12) 
with the inhomogeneous arbitrary functions 

92 U ) = iC-f), F m = l/(c?^), , F m = Gj/0%. (13) 

These space-time dependent functions obviously enter the Hamiltonian of the inhomogeneous 
quantum Toda chain as 

N 

H = Y,Pn- 9l(n + l)e U ^ U{n+1 \ Pn=Vn~ 92(n), <?l(n) = F 1{n _ 1) F 2{n) (14) 
n 

and the conserved total momentum operator as P = ^Z n P n) generated through expansion of 
the transfer matrix r(A) using the Lax operator (12) as P = CV-i, H = Cn-2 etc. 

At the classical limit one can derive the evolution equation from this inhomogeneous Hamil- 
tonian (14), which considering the relation obtained here: — ipj = iij — g^O) would recover 
the known inhomogeneous classical Toda chain equations for different choices of functions 
g a ( n ),a = 1,2 [5]. 
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4.2 Trigonometric class of inhomogeneous models: 

since both AL and SG models belong to the trigonometric class, we focus now on the idea 
behind the construction of quantum integrable models associated with the well known trigono- 
metric Rtrig- matrix [14, 19], which is a g-deformation of (2) and given through its nontrivial 
elements as 

a(A) = R{{ = Rf 2 = sin(A + 77), b(\) = R\j = R 2 2 \ = sin A, c = R l 2 \ = Rl\ = sin 77, (15) 
Following our earlier work [12] we start from a trigonometric generalization of (3): 

L ^>-{ 2sinaS+ ec 2 + e"^ 3 + r 1 ^ 53 J ' ■ ( 16 ) 

with the g-spin operators S satisfying a generalized q = e JQ -deformed algebra 



q 



S *S± = q^S±q s3 , [S\S-] = -h(cic 2 q^-c I 4q-^), h = q - q~\ (17) 



with c^a = 1,2. One can see that algebra (17) is a generalization of the well known quantum 
algebra sl q (2), which is recovered easily at the homogeneous limit cf — c 2 — 1. It is interesting 
to show that the L-operator (16) together with the _R tr j 3 -matrix (15) satisfy the quantum 
YBE (1), which yields algebra (17) as the condition for the quantum integrability. It is also 
intriguing to note that, at the undeformed limit q — > 1, the trigonometric L-operator (16) goes 
to its rational limit (3), while the quantum algebra (17) reduces to the generalized spin algebra 
(4). The construction of inhomogeneous models for this trigonometric class follows also the idea 
adopted above, i.e. we consider all central elements c^,a = 1,2, appearing in (16), (17) to be 
space-time dependent functions, which may vary arbitrarily at different lattice points j and this 
leads naturally to integrable lattice models with inhomogeneity . However since the underlying 
algebra (17) does not change even with such inhomogeneous extension, the system remains 
linked to the same trigonometric i?-matrix and hence retains its exact solvability through the 
Bethe ansatz. 

III. Q uantum inhomogeneous sine-Gordon model 

Note that an interesting realization of the quantum algebra (17) in the canonical variables: 
[u,p] = i, may be given as 

S 3 = u, S + = e~ ip g{u), S~ = g{u)e ip , (18) 

where 

g 2 (u) = — (/t + sina(s-w)(M + sina(w + s + 1) + M~cosa(w + s + 1))) (19) 
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with arbitrary constant k, spin constant s and central elements M = ±a/±1(c^C2 ± c{c 
The Lax operator (16) in realization (18) of its underlying algebra (17) would represent a 
generalized discrete sine-Gordon model, which for the trivial choice of = resulting 
M + = m 2 , M~ = 0, would recover the known exact lattice SG proposed by Korepin at al 
[15]. We introduce now the inhomogeneities by considering the central elements to be different 



8 



at different lattice points: cf = mje ia8j , cf = rrije~ iae i , which yields = m 2 M~ = 0, 
with variable mass rrij and constructs a novel inhomogeneous integrable lattice SG model linked 
to the same trigonometric R matrix [14]. It is note worthy that, in this inhomogeneous case 
we find the expression for the quantum determinant as qdetL(£) = tr(L n (^)a 2 L^(^ _1 )cr 2 ) = 
1 + Sj^q^ 1 + £~ 2 q), where Sj = (jrrijA) 2 with explicit site-dependent inhomogeneity, which 
generalizes its homogeneous case. This leads to the nonlocal interactions of the Hamiltonian 
in both the classical as well as quantum cases. Recall that in the corresponding homogeneous 
limit one gets nonlocality only in the quantum case [15]. Nevertheless, due to the property of 
the Lax-operator (16): <7 2 L*(£*)<7 2 = L(£) the Hamiltonian of the inhomogeneous lattice SG 
model as in the homogeneous case is also hermitian. 

For finding the corresponding inhomogeneous field model we may scale p and rrij by lattice 
constant A and take the limit A — > 0. This would derive from (16) through a l L t — > / + AC, 
the Lax operator C of the SG field model with variable mass parameter m = m(x,t). The 
corresponding Hamiltonian may be expressed as 

H = j dx m(x, t)(u t ) 2 + (l/m(x,t))(u x ) 2 + 8(m -m(x,t) cos(2ait)) , (20) 

describing an integrable inhomogeneous sine-Gordon field model with variable mass and placed 
in an external gauge field 6(x). Such variable mass sine-Gordon equations may arise in physical 
situations [17] and therefore the related exact results become important. 
IV. Quantum inhomogeneous Ablowitz-Ladik model 

Considering in (16) a further p-deformation of the spin operators: S ± — > = p~ s3 S ± , 
which corresponds to a p, q deformed trigonometric i?-matrix, and supposing cj~ = c\ = and 
tuning p = q , we get from the second commutator of (17) the relation 

{q^S+S- - qS~S + ) = -hc\c 2 . (21) 

If we choose c^c^ = 1 and denote S + = b, S~ = —b\ (21) reduces to a typical q-oscillator 
algebra [14]: bb' —q 2 b*b = q 2 — l, which curiously yields the known commutator in the AL model: 
[6, 6 t ] = h(l + tfb), by setting U — q 2 — 1. Note that it is easy to introduce inhomogeneity into 
the quantum AL model by considering time dependence of central elements as c\ = 4- = e tT( ^ , 

while keeping the required condition cfc~2 = 1 unchanged. This changes the Lax operator of 
the quantum AL model to its inhomogeneous version 

Lf L{mh \\) = ^ A = ^ r(t \ (22) 

which at the classical limit coincides obviously with the Lax operator of [6, 7]. Mutually 
commuting conserved operators C±j,j = N, N — 1, . . . , 1 of this inhomogeneous quantum AL 
model can be obtained from the expansion of the transfer matrix r(A), constructed using the 
Lax operator (22), while the Hamiltonian is constructed from the operators Cn-2 + C2-N 
together with a term related to the quantum determinant of (22) in the form 

H = E^-i&n + + lQg( f + h) log(l + blk), (23) 
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Here we have performed a canonical transformation b n — > b n = be l2nr W to remove explicit time 
dependence from the Hamiltonian. Though Hamiltonian (23) looks like the homogeneous one 
[14], due to b n — > b n + i2nt(t)b n , an extra inhomogeneous term would appear in the evolution 
equation and therefore using the q-oscillator type commutation relation of AL model one can 
derive from Hamiltonian (23) the known inhomogeneous discrete NLS equation at the classical 
limit. 

5. Exact Bethe ansatz solution 

Exact solutions of the eigenvalue problem for quantum integrable systems: r(A)|m >= 
A m (A)|m >, given through the algebraic Bethe ansatz (ABA) [13], can be formulated almost 
in a model-independent way. The expression for the eigenvalues may be given in the form 

A m (A) = a(A) I] M- - A) + /3(A) ft f(X - A,), where /(A) = -U (24) 

j=l j=l b K x ) 

The Bethe m-particle eigenstates are defined as \m >= YYjLi B(\j)\0 >, with the pseudovacuum 
state satisfying (7(A) |0 >= 0, where B(X) = T 12 (A), (7(A) = T 21 (A) are the off-diagonal elements 
of the monodromy matrix T(A) = Yl n L n (X). For the discrete (and also for the exact discretized 
versions of) quantum integrable models the Bethe momenta Xj are not arbitrary but should be 
determined from the Bethe equations 

(3(X 3 ) ^.a(A fe -A,)- 1 > 

Note that in (24, 25) the coefficients a (A) and /3(A) are the pseudovacuum eigenvalues of the 
diagonal elements of T(A), i.e. T n (A) | >= a(X) | >, T 22 (A) | >= /3(A) | > and 
therefore are the only model-dependent elements, since they are related to the Lax operator 
Lj(X) of a concrete model. On the other hand the factors containing the function /(A — Xj) 
, which are the major contributors in the above expressions, are given by the ratio of the R- 
matrix elements: a(A) = R\\(X) and 6(A) = R\l(X) and therefore depend not on individual 
models, but on the class to which the models belong. Consequently, they remain the same for 
all models of the same class e.g., rational, trigonometric etc. 

Therefore, we can solve the eigenvalue problem of the quantum inhomogeneous models by 
using the same formulas (24, 25) and even taking the factors /(A) same as their corresponding 
homogeneous counterparts, since in our construction we could keep the quantum i?-matrix 
same for both these cases. We have to remember however that, since the Lax operators are 
changed with the inclusion of inhomogeneity parameters, the expressions for a(X),/3(X) would 
be more complicated. Therefore for the Bethe ansatz solution of the exact lattice version of 
our inhomogeneous quantum NLS we may use (24, 25), taking the high spin limit of a(X) = 
Iln(fi , n(A + s n ) + f n ), /3(A) = Iln(fl , n(A — s n ) + f n ) , which follows from the related Lax operator 
(7). 

The model independent part for the NLS model as well as for the TC model, linked to the 
rational i?-matrix, should naturally be given by a(X) = X + a, 6(A) = A. The model-dependent 
part for the quantum inhomogeneous TC model on the other hand, as evident from its Lax 
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operator (12), is to be taken as a (A) = Y\ n (X + <72(n))> PW — 0- However the basic problem 
associated with the TC model, namely the non-availability of the pseudovacuum, remains the 
same here as in the homogeneous case. Therefore for the solution of the eigenvalue problem 
in the inhomogeneous quantum TC model one has to adopt also the functional Bethe ansatz 
method [18], in place of the above algebraic approach. 

For inhomogeneous quantum lattice SG and AL models also we can use the above Bethe 
ansatz result and since both of these models are associated with the trigonometric i? trig -matrix 
(15), we have a(A) = sin(A + a), b(\) = sin(A). For inhomogeneous SG one has to use the 
pseudovacuum by taking the product of two adjacent sites, as in its homogeneous case [19], 
though the relevant details would be different here, since the mass rrij becomes site-dependent. 

For the inhomogeneous AL model due to an additional p = g-deformation one should use 
fr±(A) = q ± b(\) in the above formulas and extra factors of q^ 1 and q should appear in the 
first and the second term of (24), respectively. Consequently another factor of q 2 also appears 
in the r.h.s. of the Bethe equations (25). The rest of the terms needed for the solution 
of this inhomogeneous AL model can be derived from its quantum Lax operator (22) giving 
a (A) = (£e ir ^) N , /3(A) = ^e tT ^)- N . 
6. Concluding remarks 

We have proposed quantum integrable extensions of the inhomogeneous NLS model, Toda 
chain and the Ablowitz-Ladik model, different forms of which are well known only at the classical 
level. We have also proposed an inhomogeneous SG model with variable mass, both at the exact 
lattice and the field limit, which are novel models at classical as well as at the quantum level. We 
have constructed such inhomogeneous quantum models exploiting the Yang-Baxter equation 
and avoiding the non-isospectrality problem by using central elements of the underlying algebra 
for introducing the inhomogeneity in a unified way. We have also indicated how to get the exact 
eigenvalue solutions of such quantum integrable inhomogeneous models systematically through 
the Bethe ansatz. It has been found that the classical limits of the quantum inhomogeneous 
models constructed here are close to known inhomogeneous equations, though their conventional 
forms can be reached only after some noncanonical transformations, permitted only at the 
classical level. The method presented here should be applicable for constructing inhomogeneous 
extension of some other quantum integrable models like the Liouville model, relativistic Toda 
chain etc. 
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